The Unitary Correlation Operator Method from a 
Similarity Renormalization Group Perspective 



O 

o 

(N 
5h . 

< 
-xt- : 

(N . 
<N : 

> : 
\o 
o 
o 
m 
o 
t> 
o 

-I— > 
I 

=5 

a 



X 



H. HergerlQ and R. RottQ 

Institut fur Kernphysik, Technische Universitdt Darmstadt, 64289 Darmstadt, Germany 

(Dated: February 4, 2008) 

We investigate how the Unitary Correlation Operator Method (UCOM), developed to explicitly 
describe the strong short-range central and tensor correlations present in the nuclear many-body 
system, relates to the Similarity Renormalization Group (SRG), a method to band-diagonalizc 
Hamiltonians by continuous unitary transformations. We demonstrate that the structure of the 
UCOM transformation, originally motivated from the physically intuitive picture of correlations in 
coordinate space, arises naturally from the SRG flow equation. Apart from formal considerations 
we show that the momentum space matrix elements of the effective interactions obtained in both 
schemes agree extremely well. 
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The key ingredient in modern ab initio studies of nu- 
clear structure are realistic nucleon-nucleon interactions, 
which reproduce the experimental two-nucleon obscrv- 
ablcs such as scattering phase-shifts with high precision. 
Different realistic interactions have been constructed us- 
ing phenomenological formulations (Argonne V18), the 
meson-exchange picture (CD Bonn), or systematic ex- 
pansions within chiral effective field theory (chiral N3LO) 
which also yield a consistent description of multi-nuclcon 
forces (for reviews, see [J Hi)- Generally, when used, 
e.g., in a no-core shell model framework these in- 
teractions exhibit very slow convergence with increasing 
model-space size The reason are large off-diagonal 
matrix elements which connect the low-energy and low- 
momentum states of the model space with high-lying 
states typically not included in the model space. In other 
words, simple model spaces do not allow for a description 
of the correlations that are induced by realistic NN in- 
teractions. 

In order to address this issue, different schemes have 
been proposed to derive soft effective interactions that 
preserve the physically constrained properties of the orig- 
inal potential, the phase-shifts. One of those, the renor- 
malization group approach [f|, introduces a momentum 
cutoff and explicitly decouples low and high-momentum 
components. Eventually one obtains a low-momentum 
interaction Ui ow -fc which is only defined below the cutoff 
scale. 

Another approach, the Unitary Correlation Operator 
Method (UCOM) 0,0, was devised as a means to de- 
scribe the dominant correlations induced by the short- 
range repulsion and the tensor interaction by an explicit 
unitary transformation. The unitary transformation of 
the Hamiltonian leads to a phase-shift equivalent cor- 
related interaction Vucom which is suitable for simple 
model spaces. It has been successfully employed in a wide 



range of many-body calculations, from the no-core shell 
model to Hartree-Fock and its extensions 0, H, 0, [l(| [HI • 
In these calculations, one benefits from the dramatically 
improved convergence behavior of the correlated interac- 
tion, which results from a prc-diagonalization of the inter- 
action in momentum space. The observed band-diagonal 
structure indicates that low and high momenta have been 
decoupled, but in a different manner than in the Uiow-fc or 
chiral EFT approach, which give just the low-momentum 
block of a block-diagonal Hamiltonian. 

Recently, Bogner et al. [H, [l]| have proposed the ap- 
plication of the Similarity Renormalization Group (SRG) 
to the NN interaction, and demonstrated the derivation 
of an effective interaction in this framework. The SRG 
approach, although starting from a different conceptional 
background, has formal similarities to the UCOM ap- 
proach. In this paper we discuss how the UCOM and the 
SRG approach can be formally related, how the resulting 
effective interactions compare, and what can be learned 
for the further development of both schemes. 

The formal structure of the UCOM is motivated by 
physical considerations on the structure and origin of the 
dominant many-body correlations. The short-range re- 
pulsion in the central part of the NN interaction drives 
the interacting nucleon pair apart. The tensor interac- 
tion induces correlations between the relative distance 
and the spin of the nucleon pair, leading to the charac- 
teristic mixing between components with relative orbital 
angular momentum L and L ± 2 in the S — 1 channel. 
To imprint these correlations on a many-body state, we 
construct a unitary transformation with the generators 



9r = (q r s(r) + s(r)q r ) 



(1) 



and 
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gn = i9(r)S 12 (r,q n ) 



#{r)-((cr 1 -qn){<T 2 -r)+ fa-r) (<r a -q n )) , (2) 



2 



where 



r 

q r - = 



2r 2 



(L x r — r x L) 



(3) 
(4) 



The generator g r uses the radial part of the relative mo- 
mentum operator q to create a shift in the radial direc- 
tion, while gtj is constructed from the so-called orbital 
momentum, i.e. the angular part of q, and generates 
shifts perpendicular to r. The strength and range of the 
transformation is governed by the central and tensor cor- 
relation functions s(r) and respectively, which are 
determined separately in each (5, T)-channel by means 
of an energy minimization 0, H| . The unitary transfor- 
mation is then written as 



C = CnC r = cxp 



cxp 



j<k 



j<k 



(5) 



where the sum runs over all nucleon pairs. One can now 
proceed to calculate expectation values either by apply- 
ing C to the many-body state or to a given observable 
O, yielding either a correlated state |^) or a correlated 
operator O: 

(*| O |#) = (tf | ClclOCnCr |$) = <tf | O |$) . (6) 

Using correlated operators is more convenient in appli- 
cations, but one has to be aware that all observables of 
interest need to be correlated consistently in that case. 

From the structure of the transformation (J5J , it is ev- 
ident that a correlated operator O will always be of A- 
body type in Fock space, where A is the number of nu- 
clcons. The decomposition into irreducible contributions 
CH"1 for a specific particle number n < A leads to the 
following cluster expansion for the correlated operator: 



o = c^oc = o [1] + o [2] 



0[A] 



(7) 



If the range of the correlation functions is small com- 
pared to the mean inter-particle distance, we can employ 
the two-body approximation and omit negligible clus- 
ter terms beyond the second order (for details see Rcfs. 
IH 0])- F° r t ne construction of the correlated Hamil- 
tonian in two-body approximation, it is then sufficient to 
consider the Hamiltonian in the two-nucleon system, 



H = T, 



int 



where we have already subtracted the center-of-mass ki- 
netic energy, which is not affected by the correlation pro- 
cedure. Applying the correlation operators, 



CtCc,HCnC T — XL t + T-' 



(9) 



and collecting the two-body contributions from the cor- 
related kinetic energy and the transformed interaction, 
we obtain the effective interaction Vucom: 



Uucom = T- 



[2] 



(10) 



Based on this definition, we can construct an explicit op- 
erator representation of Vucom , which illustrates how the 
transformation affects the structure of the interaction [H . 
The central correlator merely causes a transformation of 
the relative coordinate of the nucleon pair, 



ClrC r = R+(r), 



(11) 



where R+{r) is related to s(r) by the integral equation 



S(0 



1 . 



(12) 



The action of the tensor correlator needs to be evaluated 
via a Baker-Campbell-Hausdorff (BCH) series, which 
does not terminate in general. As a result of the Cq- 
transformation, the interaction operator acquires new 
momentum-dependent tensor terms @, which arc not 
present in typical parametrizations of realistic interac- 
tions. 

The matrix elements of Vucom in a partial-wave ba- 
sis can be evaluated directly by applying the tensor cor- 
relation operators to the two-body states, avoiding the 
complications of the BCH expansion Q. The correlated 
matrix elements in an harmonic oscillator basis have been 
used for no-core shell model calculations, which demon- 
strate the rapid convergence which results from the band- 
diagonal structure of the correlated matrix elements 

This provides a first heuristic link to the Similar- 
ity Renormalization Group (SRG) approach. Following 
Wegner's formulation of the SRG [lj|, one can write 
down a flow equation for the many-body Hamiltonian 
H . Denoting the flow parameter a, the operator evolves 



da 



[77(a), ff Q ], H Q = H, 



with 



U{a)HU\a)=T int + V a 



(13) 



(14) 



where all a-dependent contributions have been absorbed 
in the many-body interaction V a = H a — T- mt and the 
a-independent intrinsic kinetic energy T- mt — T — T cm 
has been separated. The anti-hcrmitian generator 77(a) 
formally satisfies 

r ] {a) = ^-U\a) = -r ] \a), (15) 

and has to be chosen appropriately for practical applica- 
tions. Wegner's original choice was 



77(a) = [&i&g(H a ),H a ] , 



(16) 



which one can understand intuitively: if H a commutes 
with its diagonal part with respect to a certain basis, then 
the generator vanishes and one has reached a fixed point 
of the flow. Trivial cases aside, this can only happen if 
H a is, in fact, diagonal in that basis. Thus, the generator 
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dynamically drives H a towards a diagonal structure with 
increasing a. The price one has to pay for this simplifica- 
tion is that one has to deal with complicated many-body 
interactions in V a even if one starts with a two-body po- 
tential [ilil ]. 

A simpler choice for the generator was suggested by 
Szpigel and Perry [ll| and employed successfully by 
Bogner ct al. pjj. First of all, one confines the evo- 
lution to a two-body space, discarding induced multi- 
nucleon interactions from the outset — this assumption 
corresponds to the two-body approximation used in the 
UCOM framework. The two-body generator is given by 

V(a) = [T iat ,H a } = [^,H a \, (17) 

which aims to diagonalize the two-body Hamiltonian H a 
in a basis of cigenstatcs of both p 2 and L 2 . Hence, in 
a partial- wave momentum-space basis \q(LS)JT) this 
generator drives the matrix elements towards a band- 
diagonal structure with respect to (q, q') and (L,L'). 

Let us consider the Argonne V18 [lTj potential as a 
specific example of a realistic NN interaction. Its oper- 
ator structure in a given spin-isospin (5, T) channel can 
be written as 



(18) 



where O p £ {l, L 2 , (L ■ S), (L ■ S) 2 , S 12 (£, £)}. The 
charge-dependent terms arc neglected, as they are of no 
immediate concern for the following discussion. 

Using this characteristic operator form of the interac- 
tion we can evaluate the generator (fT7|) at a = 0, which 
determines the initial flow in ([TBI: 



(19) 



v (0) = [T int ,Ho] = —[q r 
2\x 



1 ' 2 > L2 V] 



For the commutator with q 2 we find in each (S, T)- 
channcl (omitting the ST- index for brevity) , 



>^2 {q r v' p (r) + v' p (r)q r ) O p . 



(20) 



since q r commutes with all of the operators O p . For the 
angular component we get 



[E t V\ = [^ )Vt (r)SM^7)] 



(21) 



because all but the tensor component of the interaction 
commute with L 2 . Thus, we obtain the following form 
for the initial generator 



»7(0) - ^(q r S(r) + S(r)q r ) + i0(r)S 12 (r, q n ) , (22) 



with 



S(r) 



2v t (r) 
[i r 2 



(23) 

This is a remarkable result: the structure of 77(0) resem- 
bles the sum of the UCOM generators g r and gn given 
in Eqs. ([1} and @. The symmetrized radial momen- 
tum operator and the non-trivial momentum-dependent 
tensor operator S\2{r, qn), which have been constructed 
based on the structure of the short-range correlations in 
the UCOM framework, result directly from the commu- 
tation relation defining the generator of rcnormalization 
group flow. This result connects the SRC flow picture 
with a physically intuitive picture of central and tensor 
correlations in a many-body state. 

Conversely, the flow equation confirms that the UCOM 
ansatz contains the important generators. Moreover, it 
provides guidance for generalizing the UCOM scheme. 
So far, the correlation function s(r) in Eq. ([T]) describing 
the central correlations in the UCOM framework only in- 
cluded a (5, Independence. However, in the SRG frame- 
work the corresponding S(r) in (|23[) is an operator-valued 
function which is sensitive to orbital and total angular 
momentum as well. This amounts to a simple extension 
of the UCOM scheme using different correlation functions 
s(r) for each partial-wave. 

Let us return to the flow equation (fl3|) to discuss one 
inherent difference between the SRG and the UCOM 
scheme. By solving Eq. (fT3"|) , we evolve H a along a non- 
linear trajectory in the manifold of unitarily equivalent 
operators towards the fixed point, i.e. (band-) diago- 
nality. The structure of the generator 77(a) will adapt 
dynamically at each step of the flow, starting from ?;(0). 
The UCOM scheme, in contrast, uses a static generator 
to evolve the Hamiltonian in a single step from the ini- 
tial point to the final point of the flow trajectory. If one 
would assume an SRG generator independent of a, the 
flow equation could be integrated, 



(24) 



and one would recover an explicit unitary transformation 
as in the UCOM case, corresponding to a linear flow tra- 
jectory along the direction specified by 77(0). A further, 
merely technical difference is that the UCOM formula- 
tion assumes separate unitary transformations © with 
g r and go,. In the SRG picture this amounts to separate 
evolutions with two flow parameters a r and an using the 
generators 



1 1 L 2 

Vr = —[q 2 r ,Hu r ] and r)a = — [—,H ai 



2n'r 2 



(25) 



respectively. 

The difference between the static generators of the 
UCOM and the dynamic generators of the SRG affect the 
choice of the correlation functions. One could be tempted 
to use the radial dependencies of the initial SRG gener- 
ator (|2"3")) in UCOM, since they are directly given by the 
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original potential. However, this is not appropriate since 
the UCOM correlation functions effectively integrate over 
the complete flow trajectory. One would have to use the 
dynamic flow picture to determine the UCOM correlation 
functions that directly connect initial and final point of 
the flow trajectory. This is an exciting prospect and sub- 
ject of further investigations. 

The current approach for constructing the UCOM cor- 
relation functions employs an energy minimization in the 
two-nucleon system, considering the lowest partial waves 
for each(S',T)-channel[|. Let us consider the 3 Si partial 
wave with momentum eigenstates |g(01)10). Assuming 
vanishing relative momentum q, i.e. a constant radial 
wavefunction, we minimize the expectation value of the 
correlated two-body Hamiltonian 

E w = <0(01)10| ClclHCnCr |0(01)10) (26) 

by varying the parameters of suitable parametrizations of 
the central and tensor correlation functions (see Ref. Q 
for details). The range of the tensor correlation function 
is constrained using the value of the integral 



(a) AVI 8 



v(q, q') 



Iv= drr 2 §(r) . 



(27) 



The value of 1$ in the (S, T) = (1,0) channel consti- 
tutes the main parameter of Vucom, which can be un- 
derstood by scale arguments in coordinate space: the re- 
pulsive core has a very steep slope at about 0.5 fm, and 
therefore represents a distinct and physically motivated 
length scale. The case is different for the tensor interac- 
tion, which is of considerable strength up to 2 — 3 fm. 
By restricting the range of the tensor correlator, one 
introduces a scale separation — correlations at shorter 
length scales are described by the correlation operator, 
while long-range correlations need to be described by the 
Hilbert space 0]. 

That the constrained energy minimization has the de- 
sired effect becomes evident if we compare the matrix 
elements of Vucom with those of V a evolved by solving 
the full flow e quat ion (fl"3"|) in momentum space, as de- 
scribed in Ref. fl2j. In Fig. [U wc display the 3 Si matrix 
elements of AV18 along with Vucom for 1^ = 0.09 fm 3 
(optimized for calculations without 3A^-forces Q), and 



V a for a 



2 f < 



0.0215 fm 4 . This value of a was fixed 



by matching the matrix elements of Vucom an d V a at 
q = q' = Ofm -1 . 

While the matrix elements of the parent interaction 
have large off-diagonal contributions and strongly repul- 
sive diagonal terms, both Vucom and V a exhibit a dra- 
matic suppression of the off-diagonal terms and an en- 
hancement of the attractive low- momentum component. 
The Vucom and V a matrix elements are very similar, 
the only noticeable difference being the stronger increase 
of the SRG matrix elements at high momenta when ap- 
proaching the diagonal. 

As discussed above, the range constraint 1$ is the main 
free parameter of Vucom , similar to a for the SRG evo- 
lution — although this analogy should not be taken too 



q[fnr l ] 2 



(b) UCOM 
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(c) SRG 
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FIG. 1: (color online) Three-dimensional plots of the momen- 
tum space matrix elements v(q,q') in the 3 Si partial wave 
for the AV18 potential (a), the Vucom interaction with 1$ = 
0.09 fm 3 (b), and the SRG interaction V a for a = 0.0215 fm 4 
(c). The dots mark the plane of vanishing matrix elements. 



far. In Fig. [2j we compare the off-diagonal matrix ele- 
ments v(0, q) in the 3 Si and 3 Si — 3 D-y partial waves for 
various values of 1$ and a, respectively. We find that 
with increasing 1$ the decoupling of different momenta 
is improved, much like the increase of a renders V a more 
diagonal. In the 3 Si partial wave, the SRG evolution 
leads to vanishing off-diagonal matrix elements at large 
q. In comparison, the matrix elements of Vucom reveal 
tiny residual contributions at large momenta, which are 
most likely related to the particular paramctrization used 
for the UCOM correlation functions and can probably be 
remedied by using refined parametrizations or extracting 
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FIG. 2: (color online) Off-diagonal matrix elements v(0,q) 
in the 3 Si and 3 <Si — 3 D\ partial waves. Panels (a) and 
(b) show the Vucom matrix elements for 1$ = 0.03 fm 3 
( ), 0.06 fm 3 ( ), and 0.09 fm 3 ( -) in compar- 
ison to those of AV18 ( ). Panels (c) and (d) show the 

SRG-evolved matrix elements of V a for flow parameter values 

a = 0.0055 fm 4 ( ), 0.0125 fm 4 ( ), and 0.0215 fm 4 

( -), which have been chosen to reproduce the value of the 

corresponding Uucom matrix element at q = q' = Ofrn - . 

them directly from solutions of the SRG flow equation. 
Given that the correlation functions were optimized us- 
ing the diagonal 3 Si matrix element at q = q' = Ofm -1 
only, the agreement between Vucom and V a is impres- 
sive. The same observation holds for the coupled channel 
Si — 3 D±, which was not considered in the optimization 
of the UCOM correlators at all: the effects of the UCOM 
and the SRG evolution on the original interaction arc 



almost indistinguishable. Both schemes have the same 
efficiency in suppressing the off-diagonal contributions in 
this tensor-dominated partial wave. 

In summary, we have shown that the UCOM and the 
SRG, although approaching the construction of effective 
nucleon-nucleon intercations from different conceptional 
backgrounds — correlations in a coordinate space pic- 
ture for UCOM, diagonality in momentum space for the 
SRG — are closely related. Each of the methods pro- 
vide insight into the other: from the UCOM, we gain a 
physically intuitive picture for the actual effects of the 
more abstract SRG flow equation. The SRG approach, 
on the other hand, suggests ways to construct the correla- 
tion functions in a systematic fashion, without resorting 
to paramctrizations or using any parameters aside from 
the scale a, and immediately yields the structure of the 
UCOM generators. This will prove helpful when corre- 
lations in the 3/V system are studied, because the choice 
of a genuine three-nucleon correlation operator in the 
UCOM is not obvious. In the three-nucleon system, the 
SRG flow equation will likely facilitate the treatment of 
three-body correlations which are induced by the unitary 
transformation of the interaction in both approaches, and 
which have been omitted thus far. The impressive agree- 
ment of both approaches in momentum space indicates 
that the approximations required to derive the UCOM 
in its present form from the SRG — separate evolutions 
using the scale-independent generators rj r and ?yn — are 
remedied by the optimization of the correlation functions 
with physical input. 

The construction of the UCOM correlators via the 
SRG approach, as well as the construction of SRG inter- 
actions in an operator representation and the explicit in- 
clusion of three-nucleon interactions are interesting topics 
for further research. 
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